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^ ■ Abstract 

S ! The inverse medium problem for a circular cylindrical domain is studied 

^ 1 using low-frequency acoustic waves as the probe radiation. It is shown that to 

CL|| second order in k^a [ko the wavenumber in the host medium, a the radius of the 

cylinder), only the first three terms (i.e., of orders 0, -1 and +1) in the partial 
wave representation of the scattered field are non-vanishing, and the material 
\^ '. parameters enter into these terms in explicit manner. Moreover, the zeroth- 

I order term contains only two of the unknown material constants (i.e., the real 

and imaginary parts of complex compressibility of the cylinder ki) whereas the 
±1 order terms contain the other material constant (i.e., the density of the 
cylinder pi). A method, relying on the knowledge of the totality of the far- zone 
scattered field and resulting in explicit expressions for pi and ki, is devised and 
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Q ! shown to give highly-accurate estimates of these quantities even for frequencies 

^ I such that koa is as large as 0.1. 
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1 Introduction 



Retrieving the mechanical material constants (e.g., elastic moduli, density) of a material, 
either by direct measurements of these quantities, or by measurements of other variables 
from which the material constants can be derived with the help of suitable models, is one of 
the central problems in materials science. When the specimen is a solid, the elastic moduli 
are determined either by the usual static methods or by dynamic methods involving the 
inversion of data relative to resonant frequencies and/or mode shapes of vibrations excited, 
for instance by percussive forces [H], [201, [IH]; [IH]; or relative to velocities and attenuations 
for ultrasonic wave probe radiation jinj • Ultrasound methods can also be employed for fluids, 
or fluid-like materials [12], PTj . 

Another class of material characterization methods, called resonance spectroscopy j2] 
combines the underlying principles of vibratory resonances with an acoustic excitation of the 
specimen. This technique has also been employed to determine the refractive index of beads 
by means of laser irradiation pi • This class of techniques differs from the previous ones in that 
it relies on measurements of the wavefield diffracted by the specimen, and appeals to a quite 
intricate theory relating the resonances to coefficients computed from the diffracted field for 
estimating the material parameters (it can also be employed for estimating the geometrical 
parameters of the specimen jB], [2], [7j). This theory is only feasible for specimens having 
simple geometry (e.g., spherical, circular cylindrical, plate-like). A simple specimen geometry 
is also required in the standard vibration-resonance and velocity-attenuation methods if 
absolute quantitative characterizations are aimed at. 

During the last 25 years, another materials-characterization method has been developed 
which can be termed wavefield imaging. The underlying idea is: i) acquire measurements 
of the field scattered from a specimen at a series of locations in space arising either from 
several monochromatic probe fields, and/or from a pulse-like probe field, and ii) retrieve 
from these measurements an image of the specimen (i.e., a spatial map of some material 
characteristic, such as wavespeed or attenuation) with the help of a suitable model of the 
specimen/wave interaction. Insofar as there is a sharp difference between the material prop- 
erties of the specimen and those of the host medium, this method also gives a picture of 
the geometry (location, orientation, size and shape) of the specimen. When, as is often the 
case, the information relating to the material constants of the specimen is not reliable, only 
the geometrical information can been exploited (this is called qualitative wavefield imaging; 
otherwise it is called quantitative wavefield imaging). For instance, computerized diffraction 
tomography, making use of a model appealing either to the Rytov or Born approximations of 
the specimen/wave interaction, is a qualitative wavefield imaging technique except for speci- 
mens whose properties differ only slightly from those of the host medium (this is fortunately 
the case in biological imaging applications) [T7j, [121, [S], |H]- It has been suggested [T^ . 
[H], [HI that one of the reasons why Born-based techniques do not furnish reliable estimates 
of the material properties (notably the wavespeed, in specimens assumed to be lossless and 
surrounded by a host medium which is also lossless and has the same density as that of the 
specimen), is that data relating to low- frequency probe radiation was either not available or 
not used in the inversion algorithm. 

The importance of disposing of multi-frequency (and, in particular, low frequency) data 
is increasingly recognized as the key to success for material characterization in wavefield 
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imaging techniques such as the distorted Born method jTTj, |TH], the modified Born and 
modified gradient methods 0, and the contrast source method The possibihty of ob- 
taining a quantitatively-accurate image with these iterative methods is often dependent on 
being able to initialize the algorithm with a plausible image of the object at the lowest fre- 
quency of the probe radiation. More often than not, this initial image is obtained via the 
Born approximation, and since the latter is not accurate for large contrasts (between the 
host and the object) of the material constants, the algorithm has trouble with restoring the 
right values of the material constants during the iterative process. Thus, it would be useful 
to find a means for obtaining a better estimate of the material constants at low frequencies 
in the case of arbitrarily-large contrasts. This is done herein. 

In particular, we shall be concerned with the retrieval of the three material constants 
3?Ki, and pi (wherein ki = 3?ki + i'^ni) of a generally-lossy fiuid object in a lossless 
fluid host probed by plane-wave acoustic radiation. The case Q'ki = of a lossless material 
can also be treated. No assumption is made concerning the contrasts of density and com- 
pressibility between the host and the object. The latter is assumed to be a circular cylinder, 
of known radius a. The material constants of the host medium (in which the probe radiation 
propagates) are also assumed to be known, as is known the frequency and incident angle of 
the plane wave probe radiation, as well as the scattered acoustic wavefield in the far zone 
of the cylinder. The analysis for recovering the three material parameters of the cylinder is 
focused on the case in which the wavelength (Aq = ^n/ko, with ko = cu/cq, cq the speed of 
bulk waves in the host, and u the angular frequency) of the probe radiation is much larger 
than the cylinder radius. 

2 Physical configuration and governing equations 

The scattering body is an infinite cylinder whose generators are parallel to the z axis in the 
cylindrical polar coordinate system {r,6,z). The intersection of the cylinder, within which 
is located the origin O, with the xOy plane defines (see figure 1): 

i) the boundary curve T = {r = p{9);0 < 9 < 27r}, with p a continuous, single-valued 
function of 6 (further on, we shall take F to be a circle, i.e., p{6) = a), 

ii) the bounded (inner) region (i.e., that occupied by the body in its cross-section plane) 

n^ = {r < p{ey,o< e < 2tt}, 

iii) the unbounded (outer) region = {r > p{6)] < 9 < 2n}. 

It is assumed that Qq and Qi are filled with linear, homogeneous, isotropic, time-invariant 
fiuid-like media Mq and Mi respectively and that Mi is possibly lossy. 

The cylinder is probed by a monochromatic acoustic plane whose propagation vector lies 
in the xOy plane. Due to the invariance of the cylinder and incident field with respect to z, 
the incident and scattered field is also invariant with respect to z. Let U designate pressure, 
which, due to the previously-mentioned invariance, is of the form: 

f/ = f/(x,t), (1) 

with X = {x,y) = {r,6). This invariance applies also when superscripts i and d {i for 
'incident' and d for 'diffracted') are attached to U. It is convenient to associate t/(x, t) with 
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Figure 1: Problem configuration in the xOy plane 



the total field, it being understood that the latter takes the form C/j (x, t) in Q,j and: 

C/,-(x, t) = U\^, t)6jo + Uf{^, t) ; ^eQj , (2) 

with Sjk the Kronecker delta. 

We express U by the Fourier transform 

/oo 
Uj{'K,uj) exp{—iujt)duj , (3) 
-oo 

with similar expressions for W and If^. The monochromatic, plane-wave nature of the 
incident field is such that 



u 



(x, uj) = exp [-ikor cos [9 - 0')] , (4) 



wherein 9^ designates the incident angle. 

The essential task in both the forward and inverse scattering contexts is to determine 

Uj{x,u;) — u\x,uj)5jo + m^(x, a;) ; x e Qj . (5) 

via the governing equations: 

(A + A;|)ii^(x, a;) = ; X e ^- , j = 0, 1 , (6) 
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lim r^l^{dr - iA;o)M^(x, cj) = ; G [0, 27r[ , (7) 

r— >oo 

|uf(x, tu)! < cx) ; X G fii , (8) 

Mo(x, Co") — ui(r, ^, cu) = ; x G F , (9) 

aon ■ Vmo(x, uj) — aiu ■ V'Ui(x, c<j) = ; x G F , (10) 
wherein n is the unit vector normal to F and: 

k] = u\c^r\ c] = {p,K,)-\ a, = {p,)-K (11) 



3 Forward and inverse scattering problems 

The forward scattering problem (notably for simulating measured data) is formulated as 

follows: 

given: 

i) the location, shape, size and composition (material properties ) of the scattering body, 

ii) the material properties of the host medium Mq 

iii) the incident wavefield (i.e., (jH), as well as the frequency thereof, 
determine: 

the field (i.e., u'j ; j = 0, 1) scattered by the body at arbitrary points of space. 
The general inverse scattering problem is formulated as follows: 
given: 

i) the incident wavefield (i.e., (jll), as well as the frequency thereof, 

ii) the material properties of the host medium Mq 
ii) the wavefield in some subregion of f2o, 

reconstruct: 

the location, shape, size and composition of the scattering body. 

Hereafter, we shall be concerned mostly with the inverse problem, and, in particular, with 
one in which the location, size and shape of the body are known beforehand, the task being 
to reconstruct the composition of the body. 

In fact, the body will be chosen to be a homogeneous cylinder with center at the origin 
O and radius a, and we will try to reconstruct its material properties pi and ki from the 
scattered acoustic field for low-frequency probe radiation. 

4 Partial wave expressions of the fields 

A well-known identity PP informs us that the plane-wave probe radiation admits the partial 
wave expansion 

oo 

u*(x,u;) = ^ 7^ J^(fco'^) exp {im6 ) ; Vx G . (12) 
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wherein Jm{ ) is the m-th order Bessel function and 

7^ = exp {-im{d' + 7r/2) ) . (13) 
By applying the separation of variables technique to © and ((Tj), it is found that 

oo 

4{^,uj)= Cn,H^\kor)exp{tme) ; Vx G , (14) 

m=—oo 

wherein Hm{ ) is the m-th order Hankel function. The application of the same technique 
to (P) and (jHl) gives 

oo 

uf{x,uj) = DmJmikir) exp {imO ) ; Vx G f2i . (15) 

m=— oo 

From the transmission conditions Q and (fTUj) we deduce, due to the orthogonality of 
the functions {exp(im6')}, the fact that p{9) = a and n ■ V = —dr in the present case, 

Jm{koa)Jm{kia) - l3Jm{koa)Jm{kio) 

Cm = 7m ^-7Y^ TjT : , (loj 

PH^'{koa)Jm{kia) - Hfn {kQa)J„,{kia) 

Jm{koa)H^\koa) - jm{koa)H^\koa) 
(3Hin'{kQa)Jm{kia) - Hin' {kQa)Jm{kia) 
wherein: Zm{Cj = JmiO^ ^m(0 = hnear combination thereof, knowing that 

Hm\^) = JmiO + ^^^(0' with Ym{C,) the m-th order Neumann function), Z{^) := dZ{^)/d^ 
and /? = koao/kiai. 

5 Low-frequency approximation of the scattered field 
outside of the body and inversion formulas 

We use the formulas [T] 

m • 



to find: 
and 



Zo(0 = -^i(0, (19) 



-Jo{koa)Ji{kia) + (3Ji{koa)Jo{kia) 

^0 = 7o Tn m ' (20) 

-pH'i\koa)Mkia) + H^^\koa) Ji{kia) 



C±i = 

Jiikoa) [Jo{kia) - (A;ia)-Vi(A;ia)] - j3Ji{kia) [Jo(^oa) - {kQa)~^ Ji{koa)] 
71 r^T^ 77. T 77. ' (21) 



/3Ji{kia) H^^\koa) - [k^a]-^ H^^\koa) - H['\k^a) [Jo{kia) - (A;ia)-i Ji(fcoa)] 



r(i)/ 
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etc. 

We employ the notation: 

C := kilh , 5 := koa . (22) 
Due to the hypothesis of low frequencies (and/or small cylinder radius), 

5 «l . (23) 

We employ the small-argument asymptotic forms of the Bessel and Neumann functions PP 

Jo(5) ~ 1 - , Ji{5) ~ 5/2 , Yo{5) ~ (2/7r)ln(5 , Yi{5) ~ -2/715 ; 5^0, (24) 

to find 

(25) 



1715'^ 1 - C(3 



5^0 . 



(26) 



4 l + C/5 

Thus, Co and C±i are 0{5'^) as 5 ^ 0. In the same way we show that C\m\>i vanishes faster 
than 5^, so that to second order in 5 we can write 



CmH^^\kor) exp {imO ) ; Vx G ; 5^0 



(27) 



m=— 1 



or, in other terms, 

m5'^ C-/3 (1) 



4 /5 



2 1 + CP 



Hl'>{kor) cos{e~e') ; Vx G ^]o ; 5 . (28) 



It is customary, but not necessary, to measure the field in the far-field zone. In this case, we 
employ the large-argument asymptotic form of the Hankel functions PQ: 




, ^ nn TT 



^ — >• 00 



to find: 



nk^r 



exp 



TT 



kor — > 00 



(29) 



(30) 



wherein, the so-called far-field scattering function is given (in the asymptotic low-frequency 
regime) by 



iTi5'^ 



C-P ^fCP-i 



- 2 



CP+i 



cosie - 6' 



5^0 . 



(31) 



Using (HH) we find that (j2HI) and (jSH) become: 



«o(x,^) 



iT^{kQa) 



-l]H^'\kor)-2z 



a_i 

£0 

£1 + 1 

po 



H[^\kor) cos{e-e') 
; W^eQo ; (koa) , (32) 
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1l 

Ko 



- 1 




cos 



5-^0 



(33) 



The general problem in ()32|1 and ()33j) is to express A and 5 in terms of C{6) ; \f6 G [0, 27r[ 
knowing that 

C{9) = A + Bcos{9-e') ; V0G[O,27r[. (34) 



Then 



2n p2n r27Z 

C{e)cosn9de = A cosnddd + B l cos( 











9') COS n9d9 ; Vn e Z . (35) 







from which we find: 

1 Z"^'' 1 

A = — C{9)d9 , B = 

Ztt Jq tt coi 

Applied to (jH^ . this gives: 

4 1 



cos 9^ 



2n 



C{9) cos 9d9 . 



(36) 



1 + 



iTi{kQaY In Jq 



u^Q{9,9\uj)d9 



Pi = Po 



1 - . ,2 ^ rf"M;^(0,^\o;) cos^rf^ 



+ . „ ^ rf"u[^(^,^%cu)cos^d^ 

j7r(feoa) Trcost** JO UV ' ' / 



(37) 



(3^ 



This shows that: i) Ki can be retrieved independently of pi, ii) ni is a linear function of 
the measured scattered field, whereas iii) pi is a nonlinear function of this field. More 
importantly: (jH7|) - (jnH|) constitute a method for determining pi and the real and imaginary 
parts of /?! from the (far-field) scattering function, in an explicit, analytic manner. 



6 Numerical results 



We applied formulas (j?fj) and (jHHj) to retrieve the density and complex compressibility from 
the far-field scattering function. More precisely, we computed the relative errors (see figures 
2-4): 



Pi - Pi 



Pi 



(5: 



3?Ki — 3?Ki 



5. 



(39) 



(wherein pi is the value of density obtained from (jHHj) . and 3?ki, Q'ki the values of the real 
and imaginary parts of the complex compressibility obtained from (jH7|) ) over a range of 
frequencies corresponding to 10~^ < k^a < 1.5. The far-field data was simulated using 
in which the lower and upper limits of the series were replaced by -5 and +5 respectively 
and use was also made of (jHOj). The other parameters involved in the production of this data 
were: po = 1000 kg/m?, cq = 1500 m/s, pi = 1200 kg/rri^, ci = 2500 + i250 m/s, 9' = 0, 
a = 1.0 m. One sees from figures 2-4 that in order to get relative errors of the mechanical 
material parameters inferior or equal to 5%, the probe frequency should be such that koa is 
not greater than ~ 0.1. 
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Figure 2: Relative error of pi 

7 Conclusion 

Low frequency probe radiation is interesting in that it provides solutions to the inverse 
medium problem which can be written in dosed form and which are unique. Moreover, these 
solutions do not rely (as those appealing to the Born approximation) on any assumptions 
concerning the compressibility and density contrasts. Thus, the results presented herein are 
valid for arbitrary values of these contrasts. 

At the worst, the solutions of the inverse problem treated with low- frequency probe 
radiation provide suitable starting solutions for reconstructions carried out with higher- 
frequency probe radiation as well as possible explanations of the difficulties encountered in 
inverse medium problems such as the one considered herein. They may also provide decent 
estimates of the material parameters of homogeneous (and even inhomogeneous) bodies of 
more general shapes. 

In case the characteristic dimension of the body (here the radius a) is not known a priori, 
it can be determined from high-frequency probe radiation using the asymptotic technique 
described in [HI, HI [121 ■ 

The method outlined herein is obviously transposable to: homogeneous fluid slabs and 
spheres, homogeneous elastic slabs, circular cylinders and spheres, and to fluid-like or elastic 
circular tubes and spherical shells. 
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